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Abstract 

We investigate the Unruh effect for a massless scalar field in the two 
dimensional Minkowski space in the presence of a uniformly accelerated 
perfect mirror, with the trajectory of the mirror chosen in such a way 
that the mirror completely masks the Rindler horizon from the space- 
time region of interest. We find that the characteristic thermodynamical 
properties of the effect remain unchanged, i.e. the response of a uniformly 
co-accelerated Unruh detector and the distribution of the Rindler parti¬ 
cles retain their thermal form. However, since in this setup there are no 
unobserved degrees of freedom of the field, the thermal statistics of the 
Rindler particles is inconsistent with an initial pure vacuum, which leads 
us to reconsider the problem for the more physical case when the mirror is 
inertial in the past. In these conditions we find that the distribution of the 
Rindler particles is non-thermal even in the limit of infinite acceleration 
times, but an effective thermal statistics can be recovered provided that 
one restricts to the expectation values of smeared operators associated 
to finite norm Rindler states. We explain how the thermal statistics in 
our problem can be understood in analogy with that in the conventional 
version of the effect. 

PACS number: 04. 62. -|-v 


1. Introduction 


One of the central results when considering quantum field theory in the gen¬ 
eral relativistic framework is the Unruh effect, according to which a uniformly 
accelerated observer in Minkowski space perceives the Minkowski vacuum as a 
thermal bath at temperature BE] 
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where a is the proper acceleration. The effect is usually presented in two ver¬ 
sions. One refers to the distribution of particles as defined in the co-accelerated 
(Rindler) frame, which turns out to be identical with that in a thermal state of 
the field. The other considers the response of a uniformly accelerated Unruh- 
deWitt particle detector, which is also found to be of thermal form. The precise 
notion of thermality in this case is expressed by the fact that the transition rate 
of the detector respects the detailed balance or the KMS condition [3] . 

A natural question that can be asked about the Unruh effect, which led to 
various points of view in the literature, is what explains the emergence of the 
thermodynamical properties. In discussing this problem, it is important first 
to make clear what one means by “Unruh effect.” Restricting to the two ver¬ 
sions mentioned above, a useful distinction that was often made is as follows 
(although not always under these names): since the concept of particle in quan¬ 
tum field theory is a global construct, the first version can be called “global”; 
the second version, which relies on the local interaction between the detector 
and the field, can be called “local.” 

Another key concept in the discussion of the thermalization issue is that of 
the Rindler horizon, i.e. the future event horizon of the uniformly accelerated 
observer. It appears that the existing points of view can be largely distinguished 
according to the role assigned to the Rindler horizon. Simplifying the picture, 
a fundamental question to be answered is whether the thermal properties are 
intrinsically linked or not to the existence of the Rindler horizon. Without 
attempting a history of the subject, the main views along with a number of 
significant results can be summarized as follows. 

In the global version of the effect, the horizon is an essential ingredient since 
it allows to restrict the theory to the space-time region covered by the acceler¬ 
ated frame (usually a Rindler wedge), and it introduces the important distinc¬ 
tion between the observed and unobserved degrees of freedom of the field. In 
this version the thermal properties are explained in terms of the so-called ther¬ 
malization theorem unis]: the Bogolubov coefficients that mix the Minkowski 
modes and the Rindler modes have a Boltzmann form, which for observations 
restricted to a single Rindler wedge imply that the Minkowski vacuum is equiv¬ 
alent to a thermal state of the field. The mixed state implied by the thermal 
statistics is then recognized to be a consequence of ignoring the unobservable 
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degrees of freedom of the field masked by the Rindler horizon, which formally 
amounts to trace over these degrees of freedom. Thus, in this line of thought, 
thermality and the Rindler horizon are intimately linked. 

However, a different conclusion arises if one looks at the local version of 
the effect. An illuminating gedanken-experiment at this point is as follows 
[i [3 m [i uni- Consider a detector that accelerates only for a finite period 
of time, after which it follows an inertial trajectory. It is easy to see that in 
these conditions the Rindler horizon does not appear. Invoking causality and 
the locality of the field-detector interaction, one can argue that for sufficiently 
large acceleration times the transition rates will become arbitrarily close to the 
thermal rates for the infinitely accelerated trajectoriesQ This provides a simple 
example which shows that, in the local version of the effect, thermality can also 
manifest in the absence of the Rindler horizon (at least in an asymptotic sense). 

It has to be mentioned that the view that the mechanism responsible for 
the thermal rates of the detector is intrinsically local and that thermality can 
be decoupled from the existence of the Rindler horizon was firmly advocated 
more than two decades ago by Hu and collaborators in a series of papers [H 
min HU. Their conclusion was based on a calculation of the thermal rates 
using the Feynman-Vernon influence functional method which relies entirely on 
local concepts and thus emphasizes the local nature of the effect. The origin of 
thermality was explicitly identified in the exponential scaling of the quantum 
vacuum fluctuations in the detector’s proper frame, which as a purely kinematic 
effect can also manifest in the absence of the Rindler horizonjj In the author’s 
own words, “the exponential redshifting is a more basic mechanism than the 
event horizon for thermal radiance” |I4) . A somehow similar explanation in 
local terms based on the Gaussian character of the vacuum fluctuations was 
given earlier by Sciama et al. m, but without claiming the independence of 
the effect from the acceleration horizon. 

Evidence for the local nature of the effect is also provided by the recent works 
of Brenna et al. m and Brown et al. [15] , who investigated the excitation of a 

^One can explicitly see this for a variety of trajectories in m- 

^The same underlying theory was shown to apply to related phenomena such as the quan¬ 
tum radiation from accelerated mirrors and particle creation in curved space-times, with or 
without horizons; see |14II15| . 


3 



detector which accelerates inside a cavity, with different types of boundary con¬ 
ditions imposed on the field at the edges of the cavity. They found that with a 
sufficiently long interaction time and smooth switching of the interaction, the ef¬ 
fective temperature defined by the reduced density matrix of the detector shows 
a linear increase with its proper acceleration, a clear sign of the Unruh effectl^ 
The notable fact is that the effective temperature shows a negligible dependence 
on the type of boundary conditions, which again indicates a decoupling of the 
effect from the global aspects of the problem, and implicitly from the Rindler 
horizon. 

Another example which illustrates the independence of the local effect from 
the Rindler horizon was recently given by Rovelli and Smerlak uni, who con¬ 
sidered the response of a detector which uniformly accelerates in the presence 
of a perfectly reflecting plane (mirror), with the acceleration perpendicular to 
the plane. In this configuration it is possible to choose the geometry so that 
no unobservable degrees of freedom of the field exist, which makes the Rindler 
horizon irrelevant. Nevertheless, it turns out that for large acceleration times 
the transition rates become arbitrarily close to the standard thermal rates in 
empty spaceQ As stated in m, this shows that “the thermal character of the 
transition rates of the uniformly accelerated detector cannot be reduced to the 
effect of tracing out the modes behind the horizon,” which essentially coincides 
with the view in [miTaiTi] 

The aim of this paper is to add another item on the list in favour of the 
independence of the Unruh effect from the Rindler horizon. The idea in our 
construction is to mask as in m the Rindler horizon with a perfect mirror, but 
choosing this time the mirror to follow an accelerated trajectory!^ As we will 
see, a nice feature of this configuration is that it allows to analyse both the local 

^For clarity, the effective temperature obtained in EHIH] does not reproduce the Unruh 
temperature O- The distinction between the two temperatures is made clear in m- 

^An identical conclusion was previously obtained in [201121| . but without discussing the 
connection with the Rindler horizon. 

^Beginning with the paper of Deutsch and Candelas m there are a significant number of 
works which deal with quantum effects in the presence of uniformly accelerated planes (see 
e.g. |23II24| ). but these papers have little in common with our problem since they refer to the 
Rindler vacuum, whereas we consider the Minkowski vacuum. For example, a co-accelerated 
detector remains unexcited in the Rindler vacuum, which is not so in the Minkowski vacuum. 
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and global versions of the effect, much like in the usual analysis in empty space 

m- 

An aspect that somehow complicates the picture in m is that the transition 
rates depend of time and that the thermal rates are exactly obtained only in 
the limit of infinite acceleration times. It is easy to guess that the rates will be 
independent of time if the mirror follows a co-accelerated trajectory with the 
detector, in the sense that the mirror accelerates such that it appears for all 
times stationary in the detector’s proper frame. We will examine the Unruh 
effect in such a situation. For simplicity, we will focus on a massless scalar 
held in the two dimensional Minkowski space. As is well known, co-accelerated 
trajectories with uniform accelerations in the Minkowski plane can be identihed 
with hyperbolae with the asymptotes dehned by the edges of the same Rindler 
wedge. This practically hxes the kinematics in our problem. We will work as 
usual in the right Rindler wedge. A representation of the trajectories of the 
mirror and an observer/detector in the conhguration of interest is shown in 
Fig. 1. Note that for a perfect mirror all the relevant degrees of freedom of the 
held are located on the right of the mirror. It is clear that in these conditions no 
unobserved degrees of freedom exist and thus the Rindler horizon is irrelevant. 

Let us briehy state our main results. First, as expected, the transition rates 
of the detector are constant along the trajectory and differ from the rates in 
empty space. The notable fact is that they remain exactly thermal in the sense of 
the KMS condition]^ with the temperature dehned by the standard temperature 
O- Second, the distribution of the Rindler particles is also exactly thermal. 
Actually, we will hnd that the Bogolubov coefficients in the presence of the 
mirror are identical with the coefficients in empty space (for the left moving 
modes). 

®We have reported this result in a previous work m- The cited paper also contains a 
similar conclusion for a detector which accelerates in the presence of contracting and expanding 
spherical mirrors in four dimensions. 
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Figure 1: The trajectories of a uniformly accelerated mirror and of a co¬ 
accelerated observer on the right of the mirror. The dashed area represents 
the interior of the past light cone of the point P for a perfectly reflecting mir¬ 
ror. 

However, the purely thermal statistics of the Rindler particles in the presence 
of the mirror raises a problem. Since in our case there are no unobserved degrees 
of freedom of the field, the mixed state implied by the thermal statistics in 
the Rindler frame is in contradiction with a pure initial vacuum state. We 
ascribe this inconsistency to the unphysical trajectories implied by the non¬ 
zero acceleration extending to infinite past times t —>■ —oo. One should also 
recall that a well-defined initial vacuum is problematic in the presence of non- 
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stationary boundaries, which is the case here. 

Mathematically, the problem could be restated by saying that in the con¬ 
ditions of interest the formal unitary transformation which exists between the 
Minkowski and the full set of Rindler modes in empty space [5] breaks down. 
This is somehow unexpected since here too the Minkowski and the Rindler 
modes are defined on a common region of space-time in which they are both 
orthonormal and complete. We have no clear explanation for this fact. The in¬ 
consistency might have to do with the calculation of the scalar products which 
define the Bogolubov coefficients, in which we used integration by parts where 
we discarded as usual the boundary terms at infinite distances. Perhaps these 
terms are the key to the problem. In the last part of the paper we will discuss 
the thermalization of an initial well-defined vacuum, in which we will see the 
special role played by the degrees of freedom at infinite distances, which seems 
to support this idea. 

The inconsistency above invites to reconsider the problem for more physical 
trajectories of the mirror. We will naturally choose trajectories with an accel¬ 
eration that starts at some finite time to > — oo, before which the mirror moves 
inertially. We will repeat the whole analysis for these trajectories. The picture 
in the new conditions can be summarized as follows. 

Not surprisingly, the transition rates of a detector which co-moves with the 
mirror (see text) now become dependent on time, but for sufficiently large 
acceleration times they approach the thermal rates in the previous calculation. 
An analytical expression for the rates is much more difficult to obtain in this 
case, so we will be content with numerical results. As also to be expected, the 
Bogolubov coefficients for the new trajectories differ from the coefficients in the 
first calculation, depending on the time to- Given the unproblematic trajectories 
with the inertial mirror in the past, one can now expect that the Minkowski and 
the Rindler modes are related via a unitary transformation, thus preserving the 
purity of the initial vacuum in terms of Rindler states. Unfortunately, the 
complicated form of the coefficients makes a direct verification of this fact an 
extremely difficult task. 

A remarkable property of the new coefficients is that for to — oo they do 
not reduce to the old ones, which can be seen to formally correspond to the case 
to = —oo. At first sight this seems counterintuitive, but the conclusion is wel- 
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coined: if the converse were true, we would run again in contradiction with the 
purity of the initial vacuum state, which has obviously nothing to do with the 
choice of the time to- However, from a more physical point of view one can still 
expect that the thermal statistics defined by the old coefficients can be somehow 
recovered when letting to —>■ — oo. We will show that this is indeed the case. 
The key point in establishing the connection is to restrict to the expectation 
values of operators constructed from smeared creation and annihilation Rindler 
operators associated to finite norm states, which is an otherwise perfectly nat¬ 
ural restriction. This will lead to a set of “effective” coefficients which in the 
limit will reproduce the purely thermal coefficients in the first calculation. 

Finally, having in mind that in our setup no unobserved degrees of free¬ 
dom exist, the emergence of the thermal statistics from the initial pure vacuum 
requires an explanation. We will see that the restriction to smeared Rindler op¬ 
erators makes these operators insensitive to the field at large Rindler distances 
from the mirror, and we will argue that the degrees of freedom in this region 
are essential for maintaining the purity of the initial vacuum. This will lead us 
to a close analogy with the usual picture for the Unruh effect in empty space, 
in which the role of the unobserved Rindler wedge is played here by the region 
at large Rindler distances from the mirror. 

The paper is organized as follows. In the next section, we establish the form 
of the trajectories and the quantum modes in the presence of the mirror. In 
Sec. 3 we obtain the transition rates and the Bogolubov coefficients for the tra¬ 
jectories with infinite acceleration times. In Sec. 4 we obtain the corresponding 
quantities for the trajectories inertial in the past. We present our conclusions 
in Sec. 5. Technical matters are relegated to Appendices A-C. In the rest of the 
paper we use natural units such that c = Ti = ks = 


2. Preliminaries 

Let us denote by t, x the usual coordinates in the Minkowski plane. We choose 
the uniformly accelerated trajectory of the mirror to be described by 


x{t) = {f' -\- , t G (— 00 , 00 ) 


( 2 ) 



with a > 0 the proper acceleration. It is useful to introduce the Rindler coordi¬ 
nates 77 , ^ defined by 


t = — sinh ar], x = — cosh arj. 


(3) 


a a 

We recall that these coordinates cover only the right Rindler wedge x > |t|. 
The trajectory Q in Rindler coordinates is simply 


e = o. 


while the co-accelerated trajectories with dJ]) are 


(4) 


^ = constant, 


(5) 


with rj G (— 00 , 00 ). The region to the right of the mirror which is of interest 
here corresponds to ^ S (0,oo). The trajectories ([5]) in Minkowski coordinates 
as functions of the proper time r read 

t(r) = —sinh a^T, x{t) = — cosha^r, (6) 

where oj = is the proper acceleration 

The basic ingredient in our calculation are the quantum modes of the field. 
We hrst recall the form of the modes in empty space. The positive frequency 
Minkowski modes are (for both sets k G (— 00 , 00 )) 

= a; = |fc|. (7) 

VdTTWfc 

The Rindler modes are 

(^5(77,0 = c. = lkl. (8) 

VdTrwfe 

Note that the Rindler modes are defined only within the right Rindler wedge. 
Both sets of modes together with their complex conjugates are orthonormal 
and complete in their definition domain, which allows to construct the quantum 
theory with the standard procedure [I]. 

We now discuss the modes in the presence of the mirror. We ensure the 
perfect reflectivity condition by imposing Dirichlet boundary conditions (p = 0, 
which must hold everywhere along the trajectory of the mirror. The Minkowski 
modes can be readily obtained using the moving mirror model of Fulling and 
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Davies DUS]. We recall that there are two sets of inodes of special physical 
interest, i.e. the in and the out modes. From physical reasons it is clear that the 
Minkowski modes to be used in our calculation are the in modes, i.e. the modes 
which describe an unperturbed wave in the infinite past. Let us introduce the 
null coordinates 


u = t — X, V = t + X. 


(9) 


An arbitrary trajectory of the mirror can be specified via the dependence 

v=p{u). ( 10 ) 


The positive frequency Minkowski in modes for an arbitrary motion of the 
mirror in the region to the right of the mirror are D 


= uj>0. 

y/iTTUl ^ ' 

For the uniformly accelerated trajectories ([2]) the function (fTOl) 


( 11 ) 


ifl 


p(^)=-2“’ WG(-00,0). 

a^u 


We will need (ffH) in terms of the coordinates ry, Using 


a a 


( 12 ) 


(13) 


one obtains 


V dirw •' 

The Rindler modes in the presence of the mirror can be easily found since in 
Rindler coordinates the mirror appears to be stationary at ^ = 0. The positive 
frequency solutions with the desired boundary conditions are 


pSiv, 0 = -7^ , a; > 0. (15) 

V dTTo; ^ ^ 

Both sets of modes CD and CD together with their complex conjugates are 
orthonormal and complete in the region to the right of the mirror. 

We stress that the modes d and CD with r/ G (— 00 , 00 ), ^ £ (0, 00 ) are 
valid only for the trajectories (D where it is essential for t to vary over the 

^We will systematically attach a bar to the quantities defined by these trajectories. 
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entire real axis. In Sec. 4 we will deal with accelerations which begin only 
at a finite time to- The implications for the form of the modes are as follows. 
For the Minkowski modes, the new trajectories will imply a different function 
p (u) and thus a different expression (HH ( we will write the modes explicitly 
in the due place). For the Rindler modes, we can arrange things so that they 
remain unchanged. The argument is as follows. Note that in the right Rindler 
wedge the condition t > 0 is equivalent to ij > 0 and that this inequality 
selects the superior “triangular” sector of the wedge. The essential observation 
is that, if the mirror starts to accelerate at a negative time to < 0, the picture 
in the space-time sector 77 > 0 is identical with that for the trajectories with 
infinite acceleration times We will always assume that this is the case and 
concentrate on this space-time region. It is then easy to see that the same modes 
() 15 ll will appropriately describe the particles in the presence of the mirror in 
the accelerated frame. 


3. Accelerations beginning at t > — oo 

We consider in this section the case when the mirror and the detector uniformly 
accelerate for all times t G (—00, 00). We are interested in (I) the transition rates 
of the co-accelerated detector, and ( 2 ) the distribution of the Rindler particles 
in the presence of the mirror. We will use the standard Unruh-DeWitt detector 
model [I]. The co-accelerated trajectory of the detector (D) will be fixed by its 
spatial Rindler coordinate > 0 . A representation of the trajectories of the 
mirror and of the detector is shown in Fig. 2 . 
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Figure 2: The trajectories of the mirror and of a co-accelerated detector as lines 
^ = constant in Rindler coordinates. 

3.1 The response of the particle detector 

Let us denote by x(t) the trajectory of the detector in the Minkowski space, with 
T the proper time. We recall that, in the special case when the field-detector 
system is invariant to translations with respect to t, the transition rates from 
the ground state Flo = 0 to an energy level E are [T] 

/ CO 

dATe“*'®^^£)+(x(Ti),x(T2)), Ar = Ti-r2, (16) 

-OO 

where Zl''"(xi,X 2 ) is the positive frequency Wightman function in the state of 
the field. (We neglect as usual a factor that only depends on the structure of 
the detector.) The translational invariance of the system is expressed by the 
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condition 


D+{x{ti),x{t2)) = I?+(At). 


(17) 


It will become soon clear that (HB is respected in the case of interest (see eqs. 
((22l) - (l24)l below). 

As already mentioned, the role of the usual Minkowski vacuum in the empty 
space calculation is now played by the vacuum defined by the in modes (ED. 
The positive frequency Wightman function determined by these modes for an 
arbitrary trajectory of the mirror is [1] 


D^{u, v; 


, _ {p ju) - P ju') - -v' - is) 

^ ^ 'iiT ^ {v — p (u') — ie){p (m) — v ' — ie) ’ 


(18) 


For the trajectories (I12|) this translates int' 


D'^{u,v; u',v') 


= -ln(rt — u' — ie)(v — v' — ie) 

At: 


J-ln_ -- 11 ^ _. 

47r {uv'A-A/ o^ — ie){vu'+ A/a?+ ie) 


(19) 


Note that in the first logarithm we isolated the positive frequency Wightman 
function for the free field in empty space. We have to evaluate m along the 
trajectory of the detector, which is given by m where we have to set ^ 
Introducing the proper acceleration 


an = ae 


—o-Cd 


( 20 ) 


the detector’s trajectory in null coordinates is 


u{t) = - 


Qd 


;(r) = 

an 


( 21 ) 


Inserting m into m one finds 


-D'''(x(ti),x(t2)) = D^{At) +D^{At), (22) 

where (the terms and correspond to the two logarithms in (ED) 

Do (At) = In jsinh^ ~ 

®We absorbed as commonly done a positive factor into the ie term. It is essential at this 
step to use the fact that everywhere in the right Rindler wedge u < 0. 
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(At) = ^ In jsinh + a^D - ie'j sinh - a^D - ie^ | . (24) 

A convenient way to evaluate (fTHl) with the integrand defined by (l2^- (l24)) is 
as follows [3]. The first step is to eliminate the logarithms with an integration 
by parts (this can be justified with an adiabatic decoupling of the interaction at 
T —>■ ±oo), which replaces the functions In sinh z by cothx. The integrand can 
then be rewritten with 

OD ^ 

coth(^ — a) = -, (25) 

^ ^ ^ z-a + imr ^ 

n— — oo 

which allows to use a contour integration in the complex Ar plane. The inte¬ 
grals are essentially the same as in the standard calculation for the uniformly 
accelerated detector in empty space (see [25] for details). We recall that in these 
conditions the rates are (0 is the unit step function) 

K(i5) = T| (e(i5);3j4^ + e(-B)j^), = (26) 

Expressing our rates in terms of (l26l) the result is [25] 

n{E)=x{E,^D)n{E), (27) 

where 

X{E, Xd) = 1- cos {2E£,d . (28) 

Let us make a few comments on ([271). The unit term in (l28l) corresponds to 
the standard rates (IM)l . while the second term contains the perturbation due 
to the mirror. Notice the oscillatory behaviour of this term with the detector’s 
coordinate , which can be naturally seen as an interference effect between the 
incident and the reflected vacuum fluctuations in the presence of the mirror. 
(The oscillatory behaviour of the rates with respect to the detector’s position 
appears to be a constant feature in the presence of boundaries, see mmim-) 
One may also observe that the perturbation term does not vanish for infinite 
mirror-detector distances —>■ oo, which can be traced to the one-dimensional 
nature of the problem. (The perturbation terms in the four-dimensional exam¬ 
ples in [25] vanish in the same limit.) 
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It is now easy to check that the perturbed rates remain thermal. The 

thermal property in terms of the KMS condition is expressed by [3] 

n{E)/n{-E) = (29) 

with P the inverse temperature. Since (HID is respected by the thermal rates 
((^ with /3 = /3d, it is immediate from (l?7l) that the perturbed rates will respect 
the same condition provided that 

x{E,^D)=xi-E,^D). (30) 

One sees from ((Mil that relation (IM)) indeed holds. A more direct way to arrive 
to the same conclusion is to recall that the KMS condition is equivalent to [3] 
(with At seen as a complex variable) 

D+{At) = D+i-Ar-iP), (31) 


which is separately respected by (IMl) and (IMl) . 

3.2 The distribution of Rindler particles 

We now consider the distribution of the Rindler particles in the in vacuum 
in the presence of the mirror. The answer is contained in the Bogolubov coeffi¬ 
cients that mix the Rindler and the Minkowski in modes. We will focus on the 
beta coefficients, the calculation for the alpha coefficients being essentially the 
same. The beta coefficients are given by [2] 

p{I,i) = {^f*,^^), (32) 

where the parentheses denotes the usual scalar product for scalar fields. We 
first recall the result in empty space. In this case the scalar product in Rindler 
coordinates is 

f°° <->■ 

(33) 

—OO 

where the integration extends over an arbitrary hypersurface 77 = constant. Us¬ 
ing the modes in Sec. 2 and choosing rj = 0 one finds with an integration by 
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parts that for the left-moving modes k,k' < 0 the coefficients are [5] 


/ 3 (w,w') = - 


1 foJ 


27r V w' 


e 




1 

27ra 



( a \ 
W/ 




(34) 


where r(z) is the Euler Gamma function. The analogous result for the right- 
moving modes is given by the complex conjugate of (l34l) . The coefficients that 
mix the left- and right-moving modes are identically null. 

We now consider the problem in the presence of the mirror. The field in this 
case is defined only for ^ > 0 and the scalar product has to be rewritten as 


d^ip*dn‘fi2, (35) 

Jo 

where the integration can be similarly taken on an arbitrary hypersurface t] = 
constant (the vanishing of the field at ^ = 0 ensures that (|35|) remains inde¬ 
pendent of 77 ). A convenient way to perform the calculation is as follows. We 
observe that the modes in Sec. 2 can be organized as 


+ (36) 

= Gr(? + r7)-G-/(e-77), (37) 


where the form of the functions and Gui can be read from (fHl) and (fT51) . We 
further note that (the superscripts inc and ref are suppressed with no loss of 
clarity) 


d^F^{£,±'n) 

dr,GU^±Tj) 


FFfiO, 

-iuj G<^(^), 


(38) 

(39) 


where the prime denotes derivation with respect to the argument. We choose as 
before the integration hypersurface 77 = 0 in ([35]). The idea is to use integrations 
by parts in order to eliminate the derivatives introduced by (1551) . The details 
of the calculation are given in Appendix A. The result in terms of Fi^ and G^j is 

= (-2a;) x d^ {Gr(0 ^^^(0 + G™/(0 F:f/(0} • (40) 

Note that (|40l) contains no mixing between the incident and the reflected compo¬ 
nents of the modes. Using the explicit form of and G^j the integral translates 
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into 


/3(a;,a;') = - — /a) e 

2tt y uj' Jq I J 

We now change in the second term in the bracket the sign of the integration 
variable ^ —>■ —^ and observe that the two terms can be summed into a single 
integral, which leads to 

^ r~T ( /*co /aQ 


/3(w,u;') = — 

ZTT V OJ' 


1 ITT 

= - J— / (42) 

27r V w' J_oo 

Comparing with (IMl) one sees that the two integrals with respect to ^ are iden¬ 
tical. Thus the conclusion is that the two set of coefficients coincide, 


/3(w,a;') = /3(w,w'). 


(43) 


The same conclusion can be obtained for the alpha coefficients [2] 

a(/,*) = (<pf, (44) 

There is actually no need to repeat the calculation if one observes from (1321) 
and (l4^ that for the modes in the case of interest (this includes the modes in 
Sec. 4) the alpha and beta coefficients are related vig[^ 

a{uj,uj') = (45) 

It is then immediate from (1431) that (using a similar notation) 

a{ui,uj') = a{uj,uj'). (46) 


Note that thanks to (l45)) all expressions for the beta coefficients we will obtain 
here can be readily translated for the alpha coefficients. 

As a side mention, it is interesting to remark from (I40l) - (l42l) that the identity 
between the two sets of coefficients can be seen to be a consequence of the 
following facts: (i) the contribution of the incident components of the modes in 

®With the extra rule that in the factors in front of the integrals the sign of uj remains 
unchanged. 
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P{uj,uj') (i.e. the term inc-inc in (H(71) or the first integral in (H21) l is identical 
with the contribution of the region ^ > 0 in P{uj,uj'), (ii) the same is true for 
the contribution of the reflected components (i.e. the term ref-ref in (|40l) or 
the second integral in (|42])) and that of the region ^ < 0 in the same coefficients, 
and (in) there is no mixing between the incident and reflected components of 
the modes, similar to the non-mixing between the left- and right-moving modes 
in empty space. Is it as if in the presence of the mirror the “missing part” of 
the modes in the region on the left of the mirror ^ < 0 is compensated by their 
component reflected back in the region ^ > 0. 

We now recall the well-known fact mm that for the free field in the empty 
Minkowski space the form of the Bogolubov coefficients imply that the Minkowski 
vacuum appears to the Rindler observers as a thermal state with temperature 
l/P = The identities (l43l) and (I46|) show that the same is valid in the 

presence of the mirror. In particular, the distribution of the Rindler particles 
in the new conditions remains purely thermal. 

4. Accelerations beginning at to > —oo 

As we have remarked in Sec. 1, the identity between the two sets of coefficients 
raises a problem. Since in the presence of the mirror there are no unobserved 
degrees of freedom of the held, the initial Minkowski vacuum should remain pure 
when expressed in terms of Rindler states. However, this is in contradiction with 
the exact thermality implied by the coefficients in empty space. A simple way 
out of the problem is to identify the source of the inconsistency in the unphysical 
trajectories with a non-zero acceleration which extends to t —>■ —oo. We will 
thus reconsider our analysis assuming that the mirror is inertial before some 
hnite time to in the past, with to < 0 as discussed in Sec. 2. 

As a hrst step we establish the trajectories of the mirror. We assume that 
for t > to the evolution is according to the same accelerated trajectories (ED. 
We denote by J the junction point at the time t = to- We naturally choose 
the trajectory to be of class at this point. Let us denote the new trajectory 
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function (ITU)) 


V = p (u). 


(47) 


A simple calculation shows that 


p{u) = 


P-{u) if u < Uq 

p+{u) if u>uo, 


(48) 


with 



1 


(49) 


where uq < 0 is the retarded null coordinate of the point J. The connection 
with the time to is 


A diagram which illustrates the trajectory of the mirror in the new conditions 
is shown in Fig. 3 (the details on the diagram will be helpful at a later time). 
As discussed in Sec. 2 the new Minkowski in modes are given by 



(51) 


while the Rindler modes remain unchanged, 




(52) 


4.1. The response of the comoving detector 

In order to stay as close as possible to the calculation in Sec. 3.1, we will 
assume that the detector follows a comoving trajectory with the mirror. We 
mean by this that the trajectory is such that the detector appears for all times 
stationary in the proper frame of the mirror, including during the stage of in¬ 
ertial motion. Let us consider the hypersurface f} = fjo that passes through the 
point J (see Fig. 3). It is then immediate from the diagram that our condition 
implies that (1) above this line the trajectory of the detector has to coincide 

will attach a tilde to the quantities defined by the new trajectories of the mirror. 
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with the co-accelerated trajectory in Sec. 3.1, while (2) below the line the tra¬ 
jectory has to be inertial, with the velocity equal to the velocity of the mirror 
for t < to- 



Figure 3: Trajectory of the mirror with constant velocity for t < to- The lines at 
45° represent light rays. The dashed lines represent rays reflected by the mirror 
on the accelerated part of the trajectory. 

The picture above can be slightly simplified as follows. Let us denote by Tq 
the proper time of the detector at which its trajectory intersects the hypersurface 
rj = Tjo- Adopting the point of view of an observer attached to the detector, it 
is clear that in the detector’s proper frame the mirror appears to be stationary 
for T < To- Note further that by the relativistic invariance of the theory, the 
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initial vacuum in the presence of the inertial mirror coincides with the vacuum 
state in the observer’s proper frame. It is then easy to see that for t < tq the 
detector will behave like a stationary detector in the Minkowski vacuum in the 
presence of a stationary mirror (thus having a zero excitation rate), while for 
T > Tq the transition rates will depend on t only via the difference t — tq. This 
means that, without loss of generality, we can choose tq = 0. This is equivalent 
with to = 0; in which case both the mirror and the detector are fixed for all 
negative times t < 0. In this section we will focus on these trajectories. 

In the conditions above, the trajectory function (1481) simplihes to 


p{u) = 


u + 2la if u<—l/a 


— \/{a?‘u) if u>—l/a. 
The trajectory of the detector as a function of t is 

l/a_D if t < 0 


c{t) = 


+ l/a|,) 


1/2 


if t > 0, 


(53) 


(54) 


or as a function of r in null coordinates {an is defined in (1^ 1 


u{t) 


r — a 


-1 

D 


—1 f^ — CLDT 


u(r) 


T + aj^ if T < 0 

0 ) 3 ^ if r > 0. 


(55) 


It is clear that the transition rates for the new trajectories will be time- 
dependent quantities. This means that we have to consider the instantaneous 
rates. They are obtained by differentiating the transition probability along the 
trajectory with respect to t, which leads to [a [30] 

??.(£(, r) = 2 Re f dr'I?'''(x(r),x(r')) 

J —OO 

poo 

= 2Re / ds e“'®®ZI'''(x(T), x(t — s)), s = t — t' , (56) 

Jo 

where we assumed as before that the interaction begins in the infinite past. 
In the special case when the stationarity condition (na is obeyed formula (l56)l 
reduces to the time-independent rates (HSI). 

The rates in the problem of interest are in principle determined by (I56p 
together with the Wightman function pSI) . in which we have to insert the tra¬ 
jectories (1551) and ((551) . The resulting expression is rather lengthy and we do 
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not write it down here. Unfortunately, an analytical calculation seems to be 
impossible in this case, so we will be content with numerical results. Our main 
interest lies in making clear the time evolution of the rates. Also, we will only 
focus on the excitation rates of the detector, corresponding to transitions on 
positive energy levels E > 0. 

Extracting numerical values from (1561) is still not an immediate task and 
a few preliminary observations are required. An important point concerns the 
regularization of the Wightman function. It was remarked a few years ago by 
Schlicht [7] that, if one tries to obtain the rates for the uniformly accelerated 
detector using (1561) with the standard ie regularization, one does not recover 
the expected thermal rates defined by (HU). The solution proposed in [7] is to 
adopt a different regularization procedure, which essentially consists in assuming 
that the detector is an extended object with a fixed dimension in its proper 
frame. The new regularization leads indeed to the known thermal rates for 
uniform acceleratioi^ as well as to physically sensible results for various other 
trajectories HUISDlHlI The analysis in [7] considered a massless scalar field in 
four dimensions. In these conditions the Wightman function evaluated along 
an arbitrary trajectory of the detector x(t) with the new regularization method 
has the following form: 


l/47r^ 

(x(t — ie) — x(t' + ie))^ ’ 


e > 0, e —>■ 0, 


(57) 


where the squares denote the Minkowski norm x^ = and the dot 

represents derivation with respect to r. The dependence on is in the second 
fraction has to be interpreted in the sense of a formal Taylor expansion of the 
functions x(r) in powers of s, in which one retains the first two terms. 

Fortunately, there is no need to repeat the calculation in for our two- 
dimensional problem since the second formula in (1571) allows a straightforward 

Another nice feature of the new regularization is that it restores the Lorentz invariance 
of the regulated Wightman function: the fixed proper length of the detector is practically 
equivalent to a frequency cutoff in the detector’s proper frame, in contrast to the cutoff in the 
static Minkowski frame in the usual regularization. 
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generalization. The natural extension to an arbitrary Wightman function is 

£>+(x(t),x(t')) = L>+(x(t - ie),x{T' + ie)), (58) 

where in the right member £)'''(x, x') is the unregularized Wightman function. 
In the case of interest this is simply given by (fT51) with e = 0. The regularized 
form according to (1581) is obtained by making in the unregularized quantity m 
the substitutions 

(w, v) —>■ (u(t — ie), v(t — ie)), (u', v') —>■ (it(r + ie), v{t + iej), (59) 

with the same prescription for the dependence on e. The rates shown in the 
plots below are obtained from the numerical evaluation of (1561) with this regular¬ 
ization method. (The value of e was decreased until the integrals remained prac¬ 
tically unchanged; for the parameters on the plots this happens for e ^ 10“®.) 
More details on the calculation of the rates are given in Appendix B. 

The evolution of the excitation rates for a number of energies E is shown in 
Figs. 4. The vanishing values before r = 0 correspond to the interval t < 0 in 
which both the mirror and the detector are fixed. One sees that at sufficiently 
large times the rates become constant. As expected, the values at r —>• oo repro¬ 
duce the thermal rates for strictly infinite acceleration times (EZl). The curves 
in Fig. 5 represent the rates as a function of E at different times t > 0, showing 
the evolution towards the purely thermal spectrum at infinite times. The KMS 
condition for the final rates can be established from the fact that in this limit 
the Wightman function evaluated along the detector’s trajectory approaches the 
expression in Sec. 3.1, which can be easily checked with numerical calculations 
(see Appendix B). It is interesting to observe from Fig. 5 that the thermal spec¬ 
trum is roughly attained after a proper time r not very large (around ten units) 
compared to the characteristic time defined by the inverse acceleration 1/a. 
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Figure 4: The excitation rates of the detector as a function of the proper time 
T for different energies E > 0. The mirror and the detector are static up to 
the time r = 0 after which they co-accelerate as indicated in Sec. 4.1. The 
trajectory of the detector is defined by = 1. The differences between the 
constant rates at late times and the thermal rates for infinite acceleration times 
(HH) are invisible on the plot. The units for all quantities are such that a = 1 
and the same for all the following figures. 
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Figure 5: The excitation spectrum of the detector at different times r > 0 in the 
same conditions as in Fig. 4. Curve 4 (solid line) represents the purely thermal 
spectrum for infinite acceleration times (l27l) . 

Although not directly related to our subject, it is interesting to make at 
this point a comparison with the rates for the same trajectory of the detector 
in the absence of the mirror. In these conditions one has to use in (155)) the 
Wightman function for the free field, which is given by the first logarithm in 
m- It turns out that this leads to a qualitatively different time dependence of 
the rates. The essential difference is that the evolution acquires a long-lasting 
oscillatory behaviour with r. Physically, this can be seen as a consequence 
of the long range correlations for the massless field in two dimensions, which 
produce in the transition amplitude a persistent interference effect between the 
contributions due to the accelerated and the inertial part of the trajectory. In 
order to stabilize the rate, one has to eliminate the second contribution. We 
have done this by evaluating (1551) with a decoupling factor of the form e~^^, with 
77 a small positive parameter. Figure 6 shows the evolution of the rates in these 
conditions. (Compare with the times in Fig. 4.) The wild oscillations following 
the moment r = 0 are typical for a sudden variation of the kinematic parameters 
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of the detector B Note that a larger decoupling parameter ij eliminates the 
oscillations faster, but it introduces larger deviations from the expected thermal 
rates at large acceleration times. A more detailed discussion on these points is 
given in Appendix B. A notable conclusion is the mechanism which explains 
the absence of the long-lasting oscillations in the presence of the mirror: the 
Wightman function in this case is given by the difference of two logarithms (like 
in (EUl-dH), and it so happens that the oscillatory terms introduced by each 
logarithm cancel out in the difference. Intuitively, one may say that the mirror 
destroys the long-range correlations of the field along the detector’s trajectory. 



Figure 6: Evolution of the excitation rates in the absence of the mirror for two 
values of the decoupling parameter 77 (see Appendix B) and the same trajectory 
of the detector as in Fig. 4. The dashed line shows the thermal rate in empty 
space (1261) which corresponds to the uniform acceleration after r = 0. 

4.2 The distribution of Rindler particles 

^^The relevance of the low dimensionality of the problem for the oscillatory behaviour can 
be seen from a comparison with the four dimensional examples in m- For essentially the 
same trajectories of the detector the wild oscillations are absent and the evolution is very 
similar to that in Fig. 4; see the plots in the cited paper. 
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We are interested in the beta coefficients defined by the new trajectories (l48)) 

P{uj,uj') = {(p^\ (60) 

where the scalar product is identical with that in (I35p and the modes are given 
by dSB) and (l5^ . It is evident that for the Rindler modes the decomposition 
(EH) remains valid, and as can be seen from (m-® the same is true for 
the Minkowski decomposition (1361) . As a consequence we can apply the same 
calculation as in Sec. 3.2. 

We begin with the following construction. Let us denote by S the inter¬ 
section point between the hypersurface rj = 0 and the null ray reflected by the 
mirror at point J (see Fig. 3). We denote the spatial Rindler coordinate of this 
point by One can see from the diagram that the null retarded coordinate of 
S is identical to that of J, i.e. u = uq- Using (|T^ one has that 

^0 = ^ In(-aMo)- ( 61 ) 

The relation with the time tg is 

to = -^(mo + ^^o) = -- sinh(aCo)- (62) 

2 a 

Note that the limit of infinite acceleration times is equivalent to 

lim ^0 = oo. (63) 

to—^ — OO 

We now discuss the form of the Minkowski modes. Inserting (1481) in (I5ip 
and using (1131) one finds following the notation in (1551) 

Fn^ + V) = ( 64 ) 

\'4:7TCU 

and 

e-^<eo, (65) 

vdTro; 

\'4:7TCJ 

where we eliminated the parameter uq in favour of fo using m- Note that 
(l64l) and the piece of the reflected component (1651) are identical with the cor¬ 
responding quantities in the old modes ([T4)l . The difference in the new modes 
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appears only in the reflected piece (IMl) . which corresponds to the reflection on 
the inertial part of the trajectory. Also notice that the definition domains for 
(p511 and (15^ are delimited by the line ^ — r] = ^o,OT equivalently u = uq, which 
is the null ray reflected at point J. 

4.2.1 The beta coefficients 

Introducing (ITSl) and (IM1) - (I551) in (ISOl) one finds that the beta coefficients are 
(see Appendix A) 



(67) 


where 



( 68 ) 


(69) 


Before considering in more detail (I67I) - (I69I) a number of observations are neces¬ 
sary. It is quite evident that the new Bogolubov coefficients do not reproduce 
the old coefficients (1421) . Given the unproblematic nature of the new trajec¬ 
tories, one can expect that the new coefficients will preserve the purity of the 
initial vacuum in terms of Rindler states. For example, a significant difference 
between the two sets of coefficients will manifest in expectations values in the 
initial vacuum of the form ), where a^, and are the annihilation 

and creation for the Rindler particles. Since a thermal statistical matrix is di¬ 
agonal in the particle basis, the purely thermal coefficients in Sec. 3.2 imply an 
identically zero result when uji ^ 102 - As we will see in Sec. 4.3.2, this is not 
the case for the new coefficients. 

Another essential point concerns the limit of infinite acceleration times to 
—00, or equivalently ^0 —t 00. Note that the old coefficients can be re¬ 
covered from (l67l) by setting ^0 = 00 in the first integral (l68l) and ignoring the 
second integral (j69p . From the integration limits in the two integrals, one would 
expect that for ^0 —t 00 the new coefficients will reduce to the old ones. We 
will show below that this is not so. As we have remarked in Sec. I, this can be 
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seen to reflect the fact that the purity of the initial vacuum is unaffected by the 
choice of to- 

We now establish a more convenient form of the coefficients for extracting 
the limit c». Introducing the complex variables zi and Z 2 defined by 


bi : zi = — e ^ -e 

a a 


<0, 62 : Z2 = — (70) 

a a 


one can rewrite (l 68 )) and (l 6 ^ using a contour integration in the complex plan 
in the following way (we ignore a common irrelevant phase factor): 


&i(u;,cu') = --(^) “ / dti ( - - zti 

a \uj'/ Jo \a 

0 / n \ I CL 

b2(uj,uj') = X - (—-) / dt 

a \uj' J Jn 


—iu) / a— 1 


j/a—1 


0-^2 


(71) 

(72) 


Letting ^0 —t c» in the integrands the factors in the large parentheses 

vanish and one is tempted to express the integrals in terms of the Gamma 
functions T{z = =Ffa;/a). However, this is not justified due to the convergence 
condition Re z > 0, which is not respected here. The problem can be fixed by 
performing an integration by parts writing first the parentheses as a derivative, 
which increases the arguments of the Gamma functions by one. The result is 


6i(a.,a.')=--e-f« + -(4) 

LU (jJ \UJ'/ 


— iujja 


dt^ 


OJ 




— iojja 

-ih] e-‘S (73) 


62 (w, w') = H— e 

LU 


_g2*a;{o ^ 




-f-V 

LU \LU'J 


c?t 2 ( — e + it 2 


iuj j a 


r,-tl 


(74) 


The problematic limit ^0 c» reappeared in the oscillatory boundary terms 
^ but the integrals are well-defined. It is nice that in the sum (EH) 

these terms cancel each other out. 

The limit can now be easily obtained from GSl) and (EH) by letting ^0 —t 00 
under the integrals, which leads to 


lim ^(a;,a;') = {Bi{uj,uj') + B 2 {uj.,uj’)} , 

^o-s-oo ZTT V w 


(75) 


rotate the real integration axis corresponding to z\^Z 2 G [0, oo) in the original integrals 
along the negative (positive) imaginary axis in the first (second) integral. 
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where 



:“’"“/2“r(+ia;/a+1). (77) 


( 76 ) 


Note that the dependence on survives in ((77l) . Using in (176)) the recurrence 
relation 


Tiy—iL>j/a-\- 1) = {—iuj/a)T[—iijj/a)^ 


(78) 


one finds that the term Bi in (iTSl) reproduces the old coefficients (l34l) . i.e. 



(79) 


It is evident that the term B 2 does not vanish, so that it is clear that for ^0 00 

the new coefficients do not reduce to the old coefficients. 

Let us point out that in the calculation of the old/empty space coefficients, 
which are defined only by (16811 with ^0 00 , the integration by parts in (ED 

is not necessary. In this case the integral can be made convergent by simply 
adding to the exponent —ibjja — 1 a small quantity A > 0, from which ()34[) 
follows with A —>■ 0. The same procedure is inconvenient for the ^o-dependent 
coefficients, since it turns out that this mixes the dependence on the regulator 
A with that on which is undesirable. For example, repeating the integration 
by parts with an extra term A > 0 in the exponents in ED and ED one finds 
that the oscillatory boundary terms in (l7^ and (1741) do not cancel among them¬ 
selves, which introduces artificial oscillations with ^o- In addition, these terms 
acquire a factor which leads to ambiguities when trying to remove the 

regulator A —^ 0 and simultaneously let ^0 00 . 

4.3.2. Emergence of the thermal statistics for tg —00 

We have seen that for ^0 00 the term Bi in ED reproduces the empty 

space coefficients (l34l) . Hence, in order to recover the thermal statistics of the 
Rindler particles described by these coefficients it is sufficient to get rid of the 
term B 2 . The mechanism by which one can do that is rather immediate from 
the oscillatory factor in (iTTl) . The basic observation is that, in a physi¬ 

cally realistic situation, one does not measures quantities associated to a precise 
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Rindler frequency uj, but to finite norm states with a non-zero frequency width 
Aw > 0. This means that the relevant Bogolubov coefficients are obtained with 
a smearing with respect to the frequencies w. It becomes then clear that after 
integration over w the highly oscillatory factor for oo will eliminate 

the non-thermal term i? 2 - Thanks to relation (1451) the same mechanism applies 
to the alpha coefficients. This practically provides the desired connection with 
the thermal statistics obtained in the first calculation. 

To be more precise, let us denote by /n(w) the smearing functions associated 
to some set of finite norm Rindler states, with 17 standing for the mean frequency 
in the wave packet, and let us introduce the smeared annihilation and creation 
Rindler operators 

On = y dw/n(a;)a^, = J du} +. (80) 

Then the effective Bogolubov coefficients which define the measurable expecta¬ 
tion values are 

d(r7,a;') = J du) f{w) ,9(17, w') = J dw f*(u)) (81) 

Relations identical with (IMll and (IMT) can be written for the quantities in Sec. 3.2 
with the tilde replaced by the bar. According to the argument above, the 
emergence of the effective thermal statistics from the initial vacuum is expressed 
by 

lim ol(Q.,ui') = 6l{VL,ui'), lim (82) 

^ 0 —>00 

In order to illustrate the “thermalization process” implied by (15^ with some 
explicit examples, let us focus on expectation values of the form 

u(I7i,I72) = (4+h^,). (83) 

Using m together with the standard definitions of the Bogolubov coefficients 
in terms of the creation and annihilation operators mm one has that 


n(f7i, 172 ) 




(wi)/n2 (W 2 ) P{uJi,uj')P* (w2, uj') 


(84) 
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where we introduced Aa;i ,2 the frequency windows of the smearing functions. 
We recall that in empty space the thermal expectation values defined by the 
(unsmeared) beta coefficients (IMl) are [3] 

pOO 

^0 

= ( 85 ) 

As a consequence of ((43ll the same expectation values are valid for the barred 
coefficients in Sec. 3.2. The second limit in (IMll then implies 

lim n(Oi,f 22 ) 

= [ doJi f duj2 S{UJI - 0J2)fn, (c^i)/o2 (^ 2 ) —-■ (86) 

Unfortunately, due to the complicated form of our coefficients an analytical 
calculation for (1841) is virtually impossible, so we will again resort to numerical 
results. A direct numerical calculation based on (l84l) however is not simple, 
one major problem being the highly oscillatory factors when ^0 00 . 

A natural choice in order to simplify the integral is to assume very small fre¬ 
quency widths Aw, which allows to use Taylor expansions with respect to wi ^2 
in the integrand. An approximation along this line which reduces (1841) to a one¬ 
dimensional integral is presented in Appendix C. The smearing functions fn{uj) 
used in the calculation are unit-norm functions of rectangular shape centered 
in U with Aw ^ U, a (see (IC.lD b The plots below are based on the integrals 
obtained there. 

In Figures 7-10 we represented n(fli,il 2 ) as a function of ^0 for different 
values of and Aw. Note that for identical smearing functions fti = ^l 2 = ^ the 
expectation values (1841) are just the mean numbers of Rindler particles defined 
by the wave packet /n(w) in the initial vacuum, 

n(fl) = (4+h«). (87) 

This case is represented in Figs. 7 and 8. One sees that for ^0 sufficiently large 
the particle numbers approach the thermal values defined by (I86|) . i.e. 

"(») = J ■ (88) 
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The thermal numbers are not exactly recovered as a consequence of our ap¬ 
proximations (note, however, the values on the y-axis). It is interesting that 
already at = 0 (which corresponds to to = 0) the particle numbers are very 
close to the thermal values. Figures 9 and 10 illustrate the case D,i ^ 0,2- We 
have chosen the parameters in such a way that the frequency windows do not 
overlap, i.e. IIi — 172 > Aw. It is immediate from (I86p that in these conditions 

lim n(ni,r22)=0. (89) 

^ 0—>-00 

The vanishing behaviour of the expectation values when oo can be seen 

on the plots. As one could have anticipated, a larger difference fli — 172 implies 
a faster approach to zero. 





Figure 7: The number of Rindler particles in the initial vacuum n{VL) defined 
by a given fa evaluated according to Appendix C represented as a function of 
^0 • The parameters of the smearing function are shown in the box. The dotted 
line indicates the exact value defined by the effective thermal statistics (15^ . 
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Figure 8: The same as in Fig. 7 for different parameters of the smearing function. 
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Figure 9: The same as in Fig. 7 but for the expectation values n(ni, ^ 2 ) defined 
by non-identical smearing functions ^ fn^- The parameters are such that 
Q 1 —Q 2 > Aw, so that the smearing functions do not overlap. In these conditions 
the expectation values in a thermal state of the field are exactly zero. 
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Figure 10: The same as in Fig. 9 for different parameters of the smearing func¬ 
tions. 

4.3.4 An analogy with Unruh effect in empty space 

As a final point, let us present an intuitive explanation for the emergence of 
the effective thermal statistics established above. In the usual discussion of the 
Unruh effect, the fact that the vacuum appears as a mixed state for the Rindler 
observers can be understood in terms of the unobservable degrees of freedom 
of the field masked by the Rindler horizon. One can then ask what explains 
the mixed state in our case, since in the presence of the mirror no unobservable 
degrees of freedom exist. The unsmeared coefficients j5(u},u}') which do not re¬ 
duce to the purely thermal coefficients even for tg — oo can be seen to be a 
reflection of this fact. However, it is not hard to guess that the picture changes 
when we restrict to the expectation values of smeared Rindler operators. We 
will now explain that this effectively introduce unobservable degrees of freedom 
in the system, which leads to a picture very similar to that for the effect in the 
empty space. We stress that the discussion which follows is only a qualitative 
one. 
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The argument is as follows. We consider again Fig. 3. The first observation 
has to do with the state of the field on the hypersurface rj = 0, denoted by S. 
We recall that point S is the intersection of this hypersurface and the last null 
ray reflected on the inertial part of the trajectory and that the spatial Rindler 
coordinate of this point is (see (ISTjl b One can then read from the null rays 
that: (1) in the region 0 < ^ < the state of the field is locally identical with the 
perturbed vacuum in the presence of the infinitely accelerated mirror in Sec. 3.2, 
and (2) in the region < ^ < oo the state of the field is locally identical with 
the vacuum in presence of the inertial mirror. By “locally identical” we mean 
that within each region all correlation functions coincide, which follows from the 
fact that within these regions the Minkowski modes that describe the two states 
coincided Also, note that the purely thermal vacuum in Sec. 3.2 corresponds 
to the case when the local vacuum zone ^ is completely absent. This 

strongly suggest that, in order to account for the purity of the initial vacuum in 
terms of the state of the field on S, it is crucial to pay attention to the degrees 
of freedom in this region. 

The second observation is that for a smearing function of finite norm in the 
w-space the Fourier transformed wave packet in the position ^-space automati¬ 
cally vanishes for ^ oo. The important implication is that operators smeared 
with such functions will practically ignore the field at sufficiently large Rindler 
distances As a consequence, these degrees of freedom become now the “un¬ 
observed” degrees of freedom of the system. Recalling the essential role of the 
degrees of freedom at large ^ for preserving the purity of the initial vacuum 
noted above, this practically explains the emergence of the mixed state in our 
problem. The limit in which we recovered the thermal statistics was 

needed in order to ensure that the essential vacuum zone ^ is pushed suf- 
hciently far away in the region in which the smearing functions vanish for an 
arbitrary form of these functions. We are thus led to a picture very similar to 
that for the Unruh effect in empty space: the degrees of freedom of the held 
in the unobserved Rindler wedge correspond here to the degrees of freedom at 
large Rindler distances from the mirror, combined with the limit of large ac¬ 
celeration times to —^ —oo. It deserves to be remarked that the picture we 

a non-interacting theory all correlation functions are completely defined by the free 

modes. 
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have skecthed here allows to conclude that the thermal expectation values of 
the smeated operators can be also obtained for a finite tp (finite ^p), as long 
as we make sure that the support of the smearing functions in the ^-space falls 
completely outside the local vacuum region ^ > ^p. 


5. Conclusions 

Let us summarize our main results. We investigated the Unruh effect in a 
two dimensional Minkowski space in the presence of a uniformly accelerated 
perfect mirror, considering the space-time region for which the Rindler horizon 
is masked by the mirror, so that no unobserved degrees of freedom of the system 
exist. We found that the characteristic thermal properties of the effect survive, 
i.e. the response of a co-accelerated particle detector and the distribution of 
the Rindler particles in the initial vacuum still display a thermal form. Our 
result thus adds further evidence in favour of the views which emphasize the 
local nature of the effect, along with its independence from the Rindler horizon 

[6lll0l[T3[Tl]. 

For the trajectories of the mirror with accelerations which start at t > 
—oo, the transition rates of the co-accelerated detector are exactly thermal, in 
the sense that the KMS condition is obeyed. The distribution of the Rindler 
particles is also exactly thermal, being essentially identical to that in empty 
space. The last conclusion, however, is inconsistent with the absence of the 
unobserved degrees of freedom of the field in the presence of the mirror and 
the purity of an initial vacuum. We ascribed this inconsistency to the non¬ 
vanishing accelerations in the infinite past, which led us to consider the case 
when the mirror starts to accelerate at a finite time tp. 

For the new trajectories we found that the rates of a comoving detector at 
sufficiently large times approach, as expected, the thermal rates in the previous 
calculation. An interesting phenomenon emerged at a comparison with the rates 
in the absence of the mirror, in which conditions for the same trajectories of the 
detector the rates show a pronounced long-lasting oscillatory behaviour, whose 
origin can be traced to the long-range correlations for the massless field in two 
dimensions. It appears that the mirror eliminates this phenomenon. 

A notable point concerning the Bogolubov coefficients for the new trajecto- 
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ries is that for tg ^ —oo they do not reduce to the thermal coefficients in the 
previous case. We have interpreted this as a reflection of the fact that the new 
coefficients preserve the purity of the initial vacuum state in terms of Rindler 
states, irrespective of the choice of to. Nevertheless, we have shown that an ef¬ 
fective thermal statistics can be recovered for to —>■ —oo provided one restricts to 
measurements of operators constructed from smeared creation and annihilation 
Rindler operators associated to finite norm states. Formally, the mechanism is 
that the non-thermal terms in the Bogolubov coefficients contain a factor which 
in the limit becomes highly oscillatory with the Rindler frequencies, which elim¬ 
inates these terms from the smeared quantities. From a more physical point of 
view, the emergence of the thermal statistics is explained by the fact that the 
smeared operators collect information only from finite Rindler distances from 
the mirror. As a consequence of the special trajectory of the mirror, the state of 
the held within any such region for to ^ —oo becomes identical with the purely 
thermalized vacuum for inhnite acceleration times, which explains the result. 

Finally, let us mention a possible extension of our investigation. It would 
be interesting to give a more detailed description of the quantum state of the 
held on the hypersurface 77 = 0 , along with the progressive thermalization for 
to —>■ — 00 . A natural way to do this is by using the concept of entanglement 
entropy, which is generally dehned as the von Neumann entropy of the reduced 
density matrix associated to a part of the system, see e.g. [31]. In our case a 
quantity of immediate interest would be the entanglement entropy associated 
to the thermal region C < Coi or, equivalently, the entropy due to ignoring the 
vacuum zone ^ At hrst sight, one would expect to recover the Gibbs 

entropy of a thermal gas in Rindler space in the box 0 < < ^o- We want to 

point out that existing calculations in similar situations in which one restricts 
to a part of a thermal system show that the two quantities generally agree only 
in the limit of large temperatures and/or large dimensions of the box |32l 133] . 
This indicates that the entanglement entropy in our problem will reproduce the 
Gibbs entropy only in the limit of large accelerations of the mirror and/or large 
acceleration times, which probably will appear as the condition a^o ^ 1- We 
leave this as a subject for further research. 
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Appendix A 


We arrive here to the general formula for the beta coefficients (HUl) . We insert 
(I36l) - (l3^ in the integral (l3^ with rj = 0, which leads to 


poo 

Pico, w') = * / (Cr - ') 

Jo 

poo 

/ d^GT - 
Jo 


— U) X 


(A.1) 


In the first integral we eliminate the derivatives of F^^i with an integration by 
parts and replace the resulting derivatives of using (15^1) . The result is 


* r (...)= *(Gr - GIG)(F:-c ^ pref^ “ 

pOO 

-UJX / daGT-GlG){F:j^^ + Fjf). (A.2) 

Jo 

The boundary terms in (IA.2I) are zero due to the vanishing of the Rindler modes 
(the first parenthesis) on the mirror ^ = 0 and with the usual assumption that 
the field vanishes at infinite distances ^ oo. Adding the integral term in (jA.2|) 
to the second integral in (lA.ll) one obtains (HOI) . 

We now refer to the coefficients for the trajectories in Sec. 4. It is clear that 
the calculation above does not depend on the specific form of Fjj and G^j, which 
means that (HOI) can be applied to these trajectories too. The new functions 
are defined by (IM1) - (IB51) . while remain unchanged. Applying (HOI) one finds 


/3(' 


UJ.U! = —- 





V w' 1 

[^0 

Jo j 




/a) [2e-“«o_e‘*(«-2«o)] 


(A.3) 


27r V oj' 

Note that (|A.3I) differs from (1421) only via the third term which contains the 
reflected Minkowski component for 4 > ?o- Integrating in the last two terms 
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with respect to ^ —>■ — ^ and summing the first two integrals into a single integral 
one arrives to (EHi-dMl). 


Appendix B 

We detail here the evaluation of the rates (IMl)- We first make precise the for¬ 
mula used in numerical calculations. We refer throughout the appendix to the 
trajectories inertial in the past (1551) and (1551) . It is convenient to begin with 
the case of the detector in the absence of the mirror. In these conditions the 
Wightman function is given only by the first logarithm in (I19|) . which implies 
that Z)+(x(t), x(t — s)) ~ Ins for large s, so that (I5S1) is ill-defined. The con¬ 
vergence of the integral can be assured by introducing an adiabatic decoupling 
of the interaction in the infinite past via a factor e”**® with rj > 0, which allows 
to perform an integration by parts, which makes the integrand ^ 1/s. With 
these modifications the integral for the rates becomes 

= ;|lm 

^-D+(x(t), x(r - s)). (B.2) 

We ignored in (IB. Ill the divergent boundary term due to the coincidence limit 
s —>■ 0, which can be justified by the fact that for the massless field in two 
dimensions the Wightman function is defined up to an additive constant, which 
can be chosen to eliminate this term. The plots presented in Sec. 4.1 are based 
on the numerical evaluation of (jB.ip with (x, x') regularized as discussed in 
the text. For the free field and uniform accelerations the formula above with 
r] = 0 leads as it should to the standard thermal rates dMl)- 

The essential observations for the detector in the absence of the mirror are: 
(1) for large s one has W{t,s) ^ 1/s, which makes (IB. 11 ) convergent even for 
77 = 0, and (2) the integral with rj = 0 acquires a term that indefinitely oscillates 
with T. The source of this term can be understood from Fig. ll, which showJ^ 
W{t, s) as a function of s = r —r' for different values of r. The various pieces of 
the curves have the following origin: (1) the large values near s = 0 correspond 

^^Only the real part is represented. The imaginary part becomes negligible for e sufficiently 
small. The decoupling factor is not included on the plot. 


ds 


^—iEs—rjs 


W{t,s), 


(B.l) 


41 





to the coincidence limit t' —>■ t, (2) the horizontal pieces result from r' in the 
acceleration phase of the trajectory with ods ^ 1 (the Wightman function in 
this case can be approximated with (1231) where lnsinh^(oDs/2) ~ aus, and thus 
W{t,s) ~ constant), and (3) the tails at large s result from the inertial phase 
r' < 0. The oscillatory term in dHU is produced by the factor e integrated 
over the intervals which correspond to the horizontal pieces and it is directly 
related to the fact that the area below the curves indefinitely increases with r. 
The role of the decoupling factor is to eliminate this contribution. 



Figure 11: W{t,s) in the absence of the mirror represented as a function of s 
for different times r. The trajectory of the detector is the same as in Fig. 4. 
The long-lasting oscillatory behaviour of the rates in Fig. 6 results from the 
contribution of the horizontal pieces of the curves. 

A different picture arises in the presence of the mirror. The analogous curves 
are shown in Fig. 12. Note that the curves have now two peaks. The peak at 
s = 0 is the same with that in Fig. 11, while the extra peak expresses the 
large correlations between points x, x' on the detector’s trajectory which are 
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connected by a light ray reflected by the mirror. The essential fact is that 
the horizontal piece of the curves is missing, which explains why the long-term 
oscillatory behaviour does not appear in Fig. 4. The mechanism is that the 
additional logarithmic contribution besides the free field term in the Wightman 
function (1181) for large s approximately equals minus this term, which eliminates 
the horizontal piece and the tail of the curves. The rapid decrease at s —>■ oo 
leads to a well-defined rate for r —>■ oo even without a decoupling factor ij = 0. 
The rates in Fig. 6 and 7 are obtained with this value for rj. Another important 
property is that W (r, s) for r —>■ oo reproduces the quantity defined by the 
Wightman function for the uniformly accelerated trajectories (HHi-dSl , as can 
be also seen from the plot. This allows to conclude that the KMS property 
reemerges in the limit of large acceleration times. 



s 


Figure 12: The same as in Fig. 11 but for the detector in the presence of the 
mirror with trajectories of the two objects the same as in Fig. 4. The curves 
for r —>■ oo indefinitely approach curve 4 (solid line), which corresponds to the 
Wightman function for infinite acceleration times (l2^ - (l24ll . 
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Appendix C 


We arrive in this section to an approximate expression for (1841) that is well 
suited for numerical calculations. We choose the smearing functions to be unit- 
normalized rectangular wave packets of the form 


/n(w) 


l/v^ if LU G [n — rj, + rj], 

0 otherwise, 


(C.l) 


with Q the mean frequency in the wave packet and Aw = 2r] the frequency 
width. In these conditions the smeared beta coefficients are 


pQ+ri 

= / duj P{uj,uj'). 

J O —77 


(C.2) 


It is useful to parameterize the primed frequencies as uj' = ae°‘'^ with a G 
(— 00 , 00 ), which brings (15^ to the form 




dae‘^'^P{nucr)i3{il2,a)* 


(C.3) 


An obvious way to simplify (1C.21) is to assume very small frequency widths 
77 <C ill, 112 and a. Starting with (1671) an evident approximation then is (we 
replace in the front factor yCJ —>■ a/TI) 


61,2(^2, cr) 


2tt \ a 


{bi{n,a) + &2(f2,cr)} , 


(■0 + 7 ) 


Jn-r 


duj 6i,2(w, a). 


(C.4) 

(C.5) 


Similar approximations in (1C.51) must be operated with care due to the mixing 
between w and the unbounded variable ct, which can lead to non-negligible 
variations of the integrands even for small variations of w. Let us introduce 

roo 

C = a-^o, G(C,w)=/ (C.6) 

Expressing (1751) and (1751) using (|C.6I) one finds 


6l(I2, cr) 

62 ( 12 , 0 -) 


A (-O+T) 

• ^ 0+77 

xe-*-‘^G(a-eo,w), (C.8) 

v2??s2 Jn-77 
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where we eliminated the non-integral terms in (1731) and (1741) . which cancel out 
in (1C.41) . Note that we cannot make a; —>■ 14 in the factors due to the 

mixed dependence Lua. A closer look at (1C.61) shows that the same is valid in 
G((J — ^ 0 , w). In order to approximate these functions, the first step is to isolate 
the source of the rapid variations with cu. A clue is provided by the limits 


lim G(C, w) = lim G(C, w) = (c.9) 

)-oo —oo 

The expressions indicate that the rapid variations are contained only in the 
phase of G(^,a;), which is indeed confirmed by numerical calculations. With 
this simplification and considering a sufficiently small r] we can approximate 

G(C, a;)~G(C, $n(C) = ^ArgG(C, 14). (C.IO) 

Note from (jC.9p that lim^_>_oo 4)o(C) = $(14) remains finite. Introducing 
(jC.lOp in (|C.7I) and (IC.Sp the w-integrals become trivial and one finds 


6i(14, cr) 
62 ( 14 , 0 -) 



iQa „ sin[r7(o- - $n(C)] 

X -v;— 77 T-(-T 

cr - $n(C) 

^ sin[7?(2go-(j-h$n(0] 

2^0 — O' -I- $n(C) 


(C.ll) 

G(C,14). (C.12) 


It is now possible to explicitly write (jC.3p . Using (1C.41) the integral can be 
organized as (the notation is obvious) 


n — rill + ^22 + ^ 12 ) (C.13) 

where the three terms are completely dehned by (|C.11I) and (1C.121) . We write 
here only the first term: 

a r°° 

nn(14i,142)=::^^^^-^=y da G*(C, 14i) G(C, 142) 

^,^i(n,-n,)<T sin[77(a-$n,(C)] sin[7?(a - $o,(C)] 

o--$ni(C) 

The plots in Sec. 4.2.2 are based on the numerical evaluation of (jC.13l) . It is 
essential to mention that the three integrals in (1C.131) separately diverge, which 
can be seen from the fact that for a —>■ 00 the two terms (|C.7I) and (1C.81) become 
independent of a, which follows from (1C.91) . However, these divergences are just 
an artefact of our calculation. The problematic limit corresponds to w' —>■ 00 
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and going back to (1751) and (1751) one can easily check that for lo' large the beta 
coefficients behave as ^ which assures that (l84l) remains finite. In the 

numerical calculations one has to pay attention to this fact by evaluating all 
three integrals with a common upper integration limit (T^ax oo. 

We can now show that in the limit oo we recover with (|C.13I) the 

thermal expectation values defined by (1551) . Although the cr-integrals diverge, 
the fact that the total integral converges still allows to let in the integrands 
^0 oo keeping a hxed. Note that this implies —>■ —oo. It is then immediate 
from (|C.12I) that 

lim n 22 (f^i, f^ 2 ) = lim ni 2 (fli, fl 2 ) = 0. (C.I5) 

{o—>00 Co->oo 

It order to find the limit for rin, we let ( —>■ —00 in (IC.14|) and use the second 
limit in (1C.91) . which leads to 

lim fl 2 ) = , „ g-^(ni+<..)/ 2 a Y{-ini/a + l)r(i 02 /a + 1) 

7-00 cr-$(fli) cr-$(fl2) 

The trick is to recognize that the integral with respect to a is 

/ oo 1 ^00 p^li+V pf^2+ri 

da (...) = -/ da duji (C.I7) 

-00 ^ J— <x) JQi — ri 

The second relation follows from performing first the integral with respect to a 
and ignoring the second exponential factor, which is possible since the exponent 
is 00 ry. Using (1C.181) one finally finds 

lim n(IIi, 112 ) -- 1 — r(-iOi/a + l)r(ifl 2 /a + 1) 

5 o^oo 47rr7V^i^2 

rO,^+rj 

X / duji / duj 2 5 {uji — 002)- (C.19) 

J r^i —77 J ^2—^1 

It is easy to check that (1C.191) approximates the thermal expectation values 
defined by (l85l) smeared with the same functions (1C. II) . assuming a 
sufficiently small r]. The exact thermal values differ from (1C.191) only in that all 
factors above fall under the integrals with the mean frequencies ni ,2 replaced 
by the integration variables a;i^ 2 - The difference between the two quantities can 
be seen in Fig. 8 and 9. 
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